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ABSTRACT 

The ALICE experiment at LHC is dedicated to study matter formed in heavy- 
ion colhsions, but also has a strong physics program for pp coUisions. In these 
colhsions, protons wiU coUide at energies never reached before under laboratory 
conditions. At the high energies, ALICE will enable us to study jet physics 
in detail, especially the production of multiple jet events, setting the baseline 
for heavy-ion. Three-jet events allow us to examine the properties of quark 
and gluon jets, providing a suitable tool for testing QCD experimentally. We 
discuss the selection method and topology of three-jet events in ALICE. The 
analysis was performed on two PYTHIA data sets, both involving pp collisions 
at y/s — 14 TeV with enhanced jet production. The results from the dedicated 
jet MC production are discussed and compared to previous studies at CDF and 
D0. We investigate the possibilities to determine gluon jet candidates. 

1 Introduction 

In perturbative QCD, constituent partons in pp collisions experience hard scatterings. 
Outgoing scattered partons branch via soft quark and gluon radiation and this way, 
they form showers. This process is called fragmentation. The showers eventually 
hadronise and can be observed experimentally as colourless hadrons inside jet cones. 
If during the stage of fragmentation one of the partons radiates a hard gluon, one can 
expect to see a three-jet event. 



Through proper reconstruction of such events, a possibility emerges to study the 
properties of initial partons and the conditions they evolved in. Furthermore, through 
the identification of original partons, one can examine fragmentation properties of 
quark and gluon jets separately. Three jet events with the variety of processes {gg 
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999^ 11 ~^ 999^ 11 ~^ 119 those related by crossing symmetry) and topologies 
involved are suitable for thorough investigation and testing of higher order QCD. 



At LHC, protons will collide at energies never reached before and thus, opening 
doors to new physics and further scrutiny of our present knowledge. In these pro- 
ceedings we present preliminary results on topological studies of three-jet events in 
pp collisions at ^/s = 14 TeV and discuss a probabilistic approach to determine quark 
and gluon jet candidates. 



2 Quark and gluon jet differences 

Differences between quark and gluon initiated jets stem from the partons different 
colour factors, sometimes referred to as colour charges. Gluon carries colour factor 
= 3 and quark carries colour factor Cp = 4/3, hence: 

^ = - (1) 

The colour factors are proportional to the probability that a parton radiates a soft 
gluon. Colour factor of a gluon is 9/4 times bigger than the one of a quark in 
asymptotic limit of oo. This means, that gluons branch more easily and form 

higher multiplicity jets. As a consequence of the higher multiplicity, gluon jets are 
expected to be broader with softer fragmentation function than quark jets. All these 
differences can be expressed in terms of Ca/Cp ratio (Eq. [I]) [H [2]- 



Quark and gluon jet properties were previously studied by many collaborations in 
various colliding systems; ep, e^e~ and pp. Most of the studies were conducted at 
LEP e"'"e~ collider, for example see Refs. [51 S]. The source of gluon jets in these 
experiments were — > qq decays in which a hard gluon was radiated by either quark 
or anti-quark. To identify gluons in such events, the method of 6-tagging was widely 
used in two and three fold symmetry events; "Y" and "Mercedes" types respectively. 
The OPAL collaboration also performed a study with light quarks [S]. 



In e^e~ collisions, as mentioned above, the only source of gluon jets are specific 
types of events. Hadron collisions, on the other hand, involve a wider range of jet 
production processes (see Tab. [2]) and as such, provide a rich source of gluon jets. 
Studies conducted by CDF and D0 collaborations [6|,[7], inpp collisions at ^/s = 1.8 
TeV, examined quark and gluon jet properties in di-jet events. Gluon jet candidates 
were chosen on statistical basis and compared to quark jets from 7 — g and Z — q 
events. 



Experimental data collected in previous years, eg. by LEP and Tevatron collab- 
orations, indeed show different multiplicities and fragmentation properties of quark 
and gluon initiated jets. Nevertheless, as the energy scales achieved were below the 
asymptotic limit, quantitatively these differences did not reach the predictions given 
by the colour factor ratio (Eq. [1]). Study of quark and gluon jets in three jet events 
at LHC represents a further step in reaching the expected values as well as probing 
them at higher orders. 



3 Analysis 

3.1 Topology of three-jet events 

The topology of a final three parton system in hadron collisions can be described 
using eight variables. All of these are calculated in the rest frame of the final state 
three parton system. Incoming partons are labelled 1 and 2 with energies Ei and E2 
respectively (Ei > E2). Outgoing partons are labelled 3, 4 and 5 with energies E3, E4 
and E5 (E3 > E4 > E5), where parton 3 is also referred to as the leading parton. The 
partons are associated with momentum vectors p*j, i = 1,..5. Topological variables 
can be defined as follows [8], [9]: 

• Cosine of the angle between the beam and the leading parton: 

cose = —-—, (2) 

• Cosine of the angle between plane containing the leading parton 3 and the 
incoming parton 1, and the plane containing partons 4 and 5: 

, {Pi X P3) ■ {P4 X P5) 
cosip = —z ^ I I ^ ^ (3) 

\Pl X Psl-m X P5\ 

• Fractions of energies carried by each of the final state partons, called the Dalitz 
variables: 

X, = ^,t = 3A,5 (4) 

• Fractions of energies carried by a pair of partons depending on the opening 
angle between them, so called scaled masses: 



in ■ 

liij = = JXiXj (1 - cosujij) /2, (5) 



where fiij = 3,4,5) are the invariant masses of final state parton pairs, Uij 
= 3,4,5) are the angles between partons, and is the invariant mass of 
the final state three parton system. 



The two angles (Eq. [2]and[3]) determine the relative position of the final parton system 
with respect to the incoming partons (Fig. [T|, panel (b)). The energy fractions (Eq. H] 
and [SI), together with angles between partons, provide information on the event shape 
(Fig.[Il panel (a)). Experimentally, these variables are calculated for jets found by jet- 
finding algorithms. Our study is not concerned with all of the variables listed above, 
merely we focus on Dalitz variables (X3, X4) and scaled masses (/ijj, i,j = 3,4,5). 




Figure 1: Illustration of the topological variables. Panel (a): Dalitz variables and 
angles between jets, panel (b): angles with respect to incoming partons 1 and 2. 

The topology of three jet events was previously studied at Tevatron by D0 and CDF 
collaborations P, |9]. Three-jet events found in these experiments were distributed 
over the Dalitz plane with an enhancement in X3 —>■ 1 and X4 —>■ 1. With these 
values of Dalitz variables, the event shape will look like as shown in Fig. [21 in the 
[X3,X4] = [1,1] corner. Topological attributes of three-jet events in these studies 
agree with theoretical predictions. 

3.2 Sample and event selection 

Events were generated using PYTHIA 6.4 generator [10] included in the AliRoot 
framework. Samples consist of pp collisions at i/i = 14 TeV with CTEQ5L parton 
distribution functions. Events were simulated in two p^"^'"^ bins, namely p^"^^'"^ > 50 
GeV/c and p'l''"^ > 100 GeV/c, where 

phard ^j^g trausverse momentum of partons 
in the rest frame of the hard scattering pT] . 

To reconstruct jets, a UAl cone based algorithm is applied to the charged particles 
after reconstruction. More details on this jet finder can be found in [TT] and [T2]. The 
input parameters used for jet reconstruction are summarised in Tab. [H If three jets 
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Figure 2: Dalitz plane for three-jet events found using PYCELL routine inside 
PYTHIA (see text for details). In each corner of the triangle a corresponding topo- 
logical configuration is drawn. 

(each with energy > 5 GeV) are found within the fiducial region |?7| < 0.5, then 
such event is accepted. From 142000 studied events, 3080 were found to fulfil these 
conditions. 

The UAl finder does not implement splitting and merging of the found jets, there- 
fore, a particle may be assigned to more than one jet. To avoid this a separation cut 
is introduced. Separation of jets is defined as the distance between jet axes in rj — (p 
plane; AR = \J (Ar^)^ + (A0)^. All the events containing pairs of jets with separation 
smaller than twice the cone size, i.e ARij < 1.0, were filtered out. After the cut, the 
remaining sample included 1108 reconstructed events. 

The data presented are raw, without energy corrections. We compare the recon- 
structed jets to generated and PYTHIA jets. The generated ones represent jets found 
by applying UAl cone based algorithm to final state particles after generation. Same 
settings apply as for the reconstructed jets. PYTHIA jets refer to jets found by PY- 
CELL routine embedded inside PYTHIA. This routine is a cone based jet finder, that 
is applied to the final state particles after generation [TO]. The cone size used is i? = 1 
and jets within |?7| < 2.0 are accepted. 



Parameter Value 



R 


0.5 


E(seed) 


2.0 GeV 


Et ■ 

^rrnrt 


5.0 GeV 




1.5 GeV/c 



Table 1: Values of input parameters used 
for analysis. R is the cone size, E{seed) 
is the enerery of the seed, Et . is the min- 
imal energy of the reconstructed jet and 
Pf stands for minimal transverse mo- 
mentum of the tracks inside the jet. 



ID Processes ID Process 



11 


mj - 


mo 


28 


(li9 - 


(li9 


12 


QiQi - 


QkQk 


53 


99 - 


QkQk 


13 




99 


68 


99 - 


^ 99 



Table 2: 2^2 processes used in 
PYTHIA to generate QCD jets. The 
2 — i> 3 processes are generated adding 
ISR (Initial State Radiation) or FSR 
(Final State Radiation) [TO] . 



3.3 Results and discussion 

In Fig. [3] the Dalitz planes for generated (left) and reconstructed (right) three-jet 
events are shown. In the generated plot there is a clear abundance of events in the 
upper right corner of the diagram. This means that mostly events with two leading 
back-to-back jets with a perpendicular low energy jet were found by the jet finder 
(Fig. [2]). Results are consistent with CDF data. Reconstructed events show similar 




Figure 3: Dalitz plane for generated (left) and reconstructed (right) three-jet events, 
behaviour. 

The scaled mass distributions for different Dalitz variable pairs are shown in Fig. |H 
Each of the plots displays spectra for PYTHIA (open circles), generated (crosses) and 
reconstructed (full circles) events and we find good agreement between them for all 



jet pairs. Furthermore, the shapes and position of peaks are consistent with D0 data 
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Figure 4: Scaled mass distribution for pairs 34, 35, 45 in the selected three jet events. 
Open circles stand for PYTHIA jets, crosses for generated jets and full circles for 
reconstructed jets. See text for details. 



The aim of our analysis is to investigate quark and gluon jet properties and for 
that, proper determination of quark and gluon jet candidates is crucial. What we 
would like to approach is the possibility to make such a selection based on the shape 
of the event using Dalitz variables. 



The method we are going to discuss shortly is a probabilistic approach. Using this 
method we calculate the probability, that a jet carrying energy fraction Xi is a gluon, 
given other two carry energy fractions Xj, X^ (for details see Ref. [T3]). 




[1-X,){1-X^) (1-X,)(l-Xfc) (1-X,)(1- 

The probability is calculated using Dalitz variables and as such can be used to deter- 
mine gluon jet candidates based on the topology of the event (see Eq. |6]). In Fig. [5l 
probabilities of jets being gluons are shown in the Dalitz plane. We do not distin- 
guish whether a jet is quark or gluon initiated, the probability is calculated for each 
jet from the sample. Only PYTHIA jets were used. Going towards the upper right 
corner of the plane, probability that a jet is a gluon rises. This means, that this area 
of the diagram should be mostly populated by gluon jets. To prove this assumption, 
PYTHIA processes are selected from the sample. In PYTHIA, three-jet events are 
produced adding initial and final state radiation-a gluon-to the 2 —>■ 2 processes (see 
Table [2]) [10]. Most of the jets found are produced in the pure gluonic; gg gg{+g) 



(process 68), and a mixed channel; qig — > qig{+g) (process 28). This means that our 
sample is gluon dominated. Combining this fact with an abundance in the vicinity of 
the upper right corner in reconstructed Dalitz plane (contours in Fig. [5] and Fig. [3]- 
right), we can assume that the most gluons are indeed produced in area suggested by 
the probability distribution (markers in Fig. [5]) . 

pp@l4TeV, ALICE preliminary 
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Figure 5: Probability that a jet is a gluon distributed over the reconstructed Dalitz 
plane (contours) for PYTHIA jets. Empty crosses correspond to probability 0.4, 
empty rectangles to probability 0.6 and full circles to probability 0.8. See text for 
details. 




4 Summary and conclusions 

The topology of selected three-jet events and the possibility to use it in determining 
gluon jet candidates was discussed. The results on topology are presented in Figures 
[3] and HI which show Dalitz plane and scaled masses distributions, respectively. The 
shapes of the spectra are consistent with CDF and D0 results. The dominant shape 
of selected events corresponds to two back-to-back leading jets with a low energy jet 
perpendicular with respect to them (see Fig. [2]). 

To determine the gluon jet candidates, a probabilistic approach was used. Accord- 
ing to this approach, most of the gluons in the sample are produced in the vicinity of 
[X3,X4] = [1, 1] corner of the Dalitz plane. 



The discussion presented in these proceedings is based on prehminary results of the 
three-jet analysis and represents the preparation for real data analysis at ALICE. 
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